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Abstract 



Positronium annihilation amplitudes that are computed by assuming a 
factorization approximation with on-shcll intermediate leptons do not 
exhibit good analytical behaviour. Using dispersion techniques, we find 
new contributions that interfere with the known results to restore an- 
alytical properties. Those new amplitudes which cannot be obtained 
using standard factorized amplitude formalism, contribute at O (a 2 )- 
Therefore they have to be evaluated before any theoretical conclusion 
can be drawn upon the orthopositronium lifetime puzzle. 
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1 Introduction 



Positronium is a bound state of electron and positron. In this paper, we will 
be interested in the triplet state, orthopositronium, whose decay rate into 
37 has been precisely measured : 



T exp (o-Ps -► 3 7 ) = < 



' 7.0398(29) fisec" 1 Tokyo|] 
7.0514(14) fisec- 1 Ann Arbor (Gas) § 
7.0482(16) fisec- 1 Ann Arbor (Vacuum) @ 



The corresponding theoretical predictions which include perturbative 
QED corrections to a non-relativistic treatment of the bound state wave- 
function have been computed also with high accuracy (see for example 

1, 1, @): 



r (o-Ps —* 37) = am 



,.2(^- 9 ) 



.a a 2 1 / Q \2 3 a 3 1 ~ 

I -A--- \og-+B (-) -—log 2 - 
7r3 a \7r/ 2ir a 



9tt 

= (7.0382 + £ o 0.39 x 10~ 4 ) ^sec" 1 

with A = 10.286606 (10), a the fine structure constant and m the electron 
mass. Recent theoretical efforts have focused on a more complete evaluation 
of the non-logarithmic O (a 2 ) perturbative corrections, with the result B Q = 
44.52 (26) H| or, equivalently 

r {o-Ps -» 37) = 7.039934 (10) //sec -1 

This result renders the theoretical prediction still closer to the experimental 
measurement of Ref . [jl|] . 

Positronium is a test ground for bound state treatment in Quantum Field 
Theory. The first try dates back to the 40' s, with decay rates expressed 
through a factorized formula Q 

r (o-Ps -> 37) = - \(f) \ 2 ■ (Av re ia (e + e~ Sj)) v ^ 

with <j) the Schrodinger positronium wavefunction at the origin, a (e + e~ — * 3j) 
the total cross section for e + e~ — > 37 and v re i the relative velocity of e + 
and e~ in their center of mass frame. Since then, more sophisticated decay 
amplitudes have been constructed, and systematic procedures for calculat- 
ing corrections have been developed |T(J. However, the basic factorization 
of the bound state dynamics from the annihilation process has remained as 
a basic postulate. For low order corrections, this approximation is unques- 
tionable, but for O (a 2 ) corrections, factorization has to be tested. Indeed, 
non-perturbative phenomena responsible for the off-shellness of the electron 
and positron inside the positronium are of O (a 2 )- In other words, to get 
a sensible theoretical prediction at O (a 2 ), one must carefully analyze how 
binding energy effects enter the general factorization approach. 
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In a recent paper ||1 1[J , we showed that the factorization of the bound 
state dynamics from the annihilation process can be given a firm ground- 
ing through dispersion relations. From a fully relativistic model, we recov- 
ered the standard factorized amplitudes used in the literature, and found 
some forgotten O (a 2 ) contributions. Those amplitudes are gauge invari- 
ant, since they are expressed in terms of scattering amplitudes for on-shell 
e + e~ . When applied to orthopositronium, however, the appearance of those 
on-shell intermediate leptons brings infrared singularities [12|, giving an in- 
correct analytical behaviour to the amplitude o-Ps — ► 37. Indeed, basic 
principles require that this amplitude vanishes in the soft photon limit, since 
it involves only neutral bosons [13]. One must conclude that the standard 
approaches are not complete and that they can only approximately describe 
the positronium decay. 

The purpose of the present paper is to restore the analytical behaviour 
of positronium decay amplitudes. The relativistic model we used in [11] to 
describe p-Ps — ► 77 has the correct analytical behaviour when applied to 
o-Ps — > 37 since only off-shell constituents appear. Indeed, we will find 
along with the known factorized-type amplitudes already found for para- 
positronium, a whole class of non-factorizable processes which enforce the 
positronium decay amplitude to vanish in the soft photon limit. Those new 
amplitudes were absent for p-Ps — ► 77, but will contribute to o-Ps — ► 37 at 
0{a 2 ). 

This paper is organized as follows. We first give a brief description of 
the derivation of factorized amplitudes from our relativistic model. Then we 
analyze the simple decay p-Dm — ► je + e~ , and show in detail how the ana- 
lytical behaviour of the amplitude is restored by new contributions. Finally, 
we discuss the orthopositronium decay to three photons and describe all the 
new contributions. In the appendix, we present a short discussion of the de- 
cay Kg — ► 7e + e~, the hadronic counterpart of the decay p-Dm — > "ye + e~ , in 
order to show the generality of the arguments developed for electromagnetic 
bound states. 



2 Factorized Amplitude from a Loop Model 

Here we briefly review the recently proposed method [jll| to derive factorized 
amplitudes for bound state decays. 

The decay rates for positronium are calculated in a loop model. The 
positronium first decays into a virtual electron-positron pair, which subse- 
quently annihilates into real or virtual photons (an odd number for ortho- 
states, an even number for para-states). The coupling of the positronium 
to its constituents is described by a form factor, denoted by Fb- It is not 
assumed to be a constant, since a constant form factor would amount to 
consider positronium as a point-like state. 
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In our model, the decay amplitudes are (see for example figure 4) 

(2) 

for parapositronium (pseudoscalar) and orthopositronium (vector, with po- 
larization e a ) respectively . The T^ v "' is the scattering amplitude for off- 
shell e + e~ with incoming momenta ^P + q and \P — q into the final state 
A, m is the electron mass, P the positronium four- momentum and Fb = 
Fb (q 2 ,q • P). Prom these amplitudes, the decay widths are calculated as 

r (^S - A) = — - — [ d<S> A Y \M^- ( 2J+1 S - A) e*el...\ 2 

pol 

with M the positronium mass. 

What we have shown in our previous article [|ll| is that when one ex- 
presses the loop integration in ([l]) or (|2|) via a dispersion relation (fli|], 
fll5| ), one ends up with a factorized convolution- type amplitude, i.e. for 
parapositronium: 

M »u... { i s ^ A) = _J ( k 2) Tr | 7s (-}6> + m )r^-{k,k',h,...)Qt + m 

(3) 

with if) (k 2 ) the bound state wavefunction, r^"*" (k, k', l±, ...) the reduced 
scattering amplitude for on-shell e + e~ with momenta k and k', k = — k' 
and C = \[M jm (M is the positronium mass). The analogous form for 
orthopositronium decay is obtained using the substitution 7 s — >0. The form 
factor Fb is related with ip (k 2 ) through 

F B = C^(k 2 ) (k 2 + 7 2 ) 

with O ^" (k 2 ) = ip (k 2 ), O being the Schrodinger wavefunction at the ori- 
gin, and 7 2 = m 2 — M 2 /4 (related to the binding energy and fine-structure 
constant through Eb = M — 2m = —ma 2 / 4). 



3 Paradimuonium Decay into ^e + e 

The present section concerns paradimuonium, the singlet fjr fjT bound state 



[ 16 1 . This state has not been observed yet. The reason to consider the decay 
p-Dm — > 7e + e~ is that it is the simplest process where a photon is kine- 
matically allowed to have a vanishing energy. Therefore, we will be able to 
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test whether the proposed factorization procedure (Q) gives a correct ana- 
lytical behaviour to the amplitude in the soft photon limit. The simplicity 
of the process comes from the pseudoscalar nature of the paradimuonium, 
which allows a manifestly gauge invariant treatment throughout. Note that 
the positronium decay p-Ps — ► "yvV via a Z° has the same dynamics. Af- 
ter having analyzed this simple decay, we will be ready to tackle the more 
interesting process o-Ps — > 777. 

The decay p-Dm — > e + e~7 is shown in figure 1. When the form factor 
Fb appearing at the vertex p-Dm — > p + p~ allows a change of variable 
q — ► —5, the two amplitudes can be combined into 

M {p-Dm - e + e- 7 ) = SmeV^V^* (jfe) i^Ii^Mlj ( M 2 ,x) 

(4) 

with x = 2to /M the reduced photon energy, M the dimuonium mass and m 
the muon mass. The loop integral, which can be viewed as an effective form 
factor, is given by 

Jfp2 x) / 1 1 1 

^ ' J J (2vr) 4 ( g -ip) 2 - m 2 (g+ip) 2 -m2 (g-iF + fc) 2 -r 

where 77 = P 2 /M 2 and f 2 = P 2 (l — x). Gauge invariance is ensured by the 
factorized tensor structure, i.e. by the antisymmetric Levi-Civita tensor. 
I (P 2 ,x) is given by the same expression as the two-photon decay integral 
we encountered previously [11], the only difference being that t 2 7^ 0. 
The differential width is given by 

dT (p-Dm — > e + e~7) 16a 3 o, r r>,, r n N ,2 , . 

^- - = m 2 M 3 \l(M 2 ,x)\ p(x,a) 

dx 3 

with the phase space factor 



p(x,a) = *ll [a + 2 (1 — x" 
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2 



1-z (1-z) 

where a = 4m 2 /M 2 , and the bounds on x are [0, 1 — a]. 

We are going to calculate the integral X(M 2 ,x) using dispersion rela- 
tions. The main difference with the two-photon decay case is the appearance 



of two different cuts (figure 2). We have shown [11| that considering the ver- 
tical cuts is strictly equivalent to the decay amplitude calculation done using 
formula (@) where T^ is the scattering amplitude p + (fe') /j - (k) — > 77* — > 
e + e _ 7 with on-shell muons. Standard approaches found in the literature are 
then obtained by making an O (a 2 ) approximation in (^). Obviously, those 
approaches completely miss the oblique cuts. What we are going to show is 
that those forgotten cuts contribute at the order of 7 2 = m 2 — M 2 /A ~ a 2 /4, 
and more importantly, that the amplitude for p-Dm — ► e + e~7 has a correct 
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soft-photon limit behaviour only if we take all the cuts into account. By this 
we mean that each cut gives a contribution to the amplitude that behaves 
as a constant when the photon energy goes to zero. The combination of 
the vertical and oblique cuts, on the contrary, forces the amplitude to van- 
ish in that limit. We thus recover the analytical behaviour expected from 
Low's theorem [O] for the decay of p-Dm — ► 77* involving only neutral 
bosons (for p-Dm —* je + e~ , infrared divergences cannot arise from electron 
bremsstrahlung processes due to selection rules). 



3.1 Soft photon limit property of the absorptive part 

To keep the discussion as general as possible, we extract the imaginary part 
of I (P 2 ,i 2 ) without specifying the form factor Fg. The absorptive part is 
found by cutting the relevant propagators as 



2 2 
— m 



2(2^) 4 ( g _i P + ^_ m 2 

4 2vrz(5 (" (q + \P) 2 - m 2 ) 2vr^ ( (g - |P + jfe) 



2(2*0* (q-\P)' 



for the vertical and oblique cuts, respectively. By a straightforward integra- 
tion, the first expression gives (s = P 2 ): 



„Fb (go = 0, |q| = yV 4 - m2 
ImXi (s, x) = — In 

S 16lTX 



1 + ^1 - Am 2 /s 
1 - y 7 ! - 4m 2 /s 



fs - 4m 2N ) 



(5) 

while the second one cannot be completely integrated without specifying 
F B - 



Tj 1 



<?max F B (g , |q| = s/q' 2 + qo^/s + s/4 - m 2 ) / 4m 2 



ImJ 2 (s,x) = -— — / (fg -0 s 

s 16?rx J qui . n q \ 1-x 

(6) 

with the bounds given by 



X\fs I / Am 2 1 \ X\fs I / 4m 2 1 

<?min = T~ 1 + V 1 1 i <?max = r" 1 ~ V 1 



8 1 — X I 4 \ V 5 I — X 

(7) ' 

Interestingly, the second cuts contribute for go 7^ 0- This is in sharp con- 
trast with the two-photon decay, since there only the first cuts exist. Fur- 
ther, approximated bound state wavefunction where a J (go) appears cannot 
be used ([5], ||, 0), and one should revert to the full four-dimensional 
Bethe-Salpeter wavefunction (for example the Barbieri-Remiddi one [17]). 
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Let us now demonstrate an important point, i.e. that in the soft photon 
limit, the combination ImXi (s, x) + ImX2 (s, x) behaves as a constant when 
x — > despite the fact that each cut diverges in that limit. This will 
guarantee that the imaginary part of the whole amplitude vanishes in the 
soft photon limit thanks to the presence of k p in the tensor structure of (|j). 
This behaviour is what one expects from Low's theorem For all the 

form factors we considered, the behaviour of the combination ImXi (s,x) + 
ImX2 (s,x) when x — > is summarized by : 



ImXi (s, x) + ImX 2 (s, x) --7— 

S lOTT 



F B (0) | Vs - 4m 2 8F B 



1 



4m 2 
s 



dq 



Specifically, the form factors of interest are constructed from the Schrodinger 
wavefunction (see [0]). Let us recall their definition and give the limit 
obtained when using each of them : 



Fb = Wo 



8-7T7 

q 2 + 7 2 



s 2 



(8) 



ImXi (s, x) + I111I2 (s, x) x — >° — ^° 7- 



m 



1 



4m 



7 Z 



m 2 + - 



Fb 1 = C0 O 



32^7^ 
(q 2 + 7 2 ) 2 



0) 



ImX, (,,*) + ImX 2 (s,x) ^ (V ~ ' ( * ~ 



s 2 



1 



4m 



7 2 — m 2 H — 



i.e. constant limits. Remark that for a constant form factor Fb = Xg onst , 
the limit is simply 



ImXi (s, x) + ImX2 (s, x) 



7711 



Const 



1 



167T y/l - 4m 2 /s 



It is easy to verify that this result is indeed what can be calculated from the 

imaginary parts : 

ImXi (s, x) + ImX2 (s, x) = 



V F i 



Const 
B 



S WlTX 



hi 



i-a/i- 



4m^ 



(s — 4m 2 ) - In 



s(l-a:) 



l-l/l- 



4m^ 
s(l-3;) 



In conclusion, the absorptive part of the amplitude is seen to vanish. 
Thus, we expect that the dispersive part will also vanish in that limit. The 
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oblique cuts, omitted from standard analyzes, are seen to be essential to 
maintain good analytical properties of the amplitude. 

In the appendix, the restoration of analyticity by a cancellation among 
the vertical and oblique cuts is presented in the context of the kaon decay 
K% — ► e + e~7 via a charged pion loop. In this analysis, the form factor 
describing the K — ► tttt vertex is taken as a constant. This decay is inter- 
esting since there is a close similarity between this hadronic decay process 
and the present QED bound state decay. The main conclusion is that even 
if charged particles are present in intermediate states, the amplitude has to 
vanish in the soft photon limit, as expected from the fact that the decays 
K s — ► 77* or p-Dm —* 77* involve only neutral bosons. Further, this K s 
decay provides a physically sensible process where one can analyze the con- 
stant form factor assumption, since for such a loosely bound system as the 
p-Dm a constant form factor cannot be realistic. 

We will now evaluate the contribution of each cut for the case of the 
Schrodinger momentum wavefunction form factor in (|| ). In doing so, 
we will gain a better understanding of the analyticity restoration at the 
level of the differential decay rate for p-Dm — ► e + e~7. The purpose of the 
next sections is only illustrative, and the discussion can be straightforwardly 
transcribed for the advocated improved form factor F^ 1 (g) (see [11]). 



3.2 The vertical cuts reproduce standard approach results 

To precisely show what happens when the oblique cuts are forgotten, we 
compute here the rate with the vertical cuts only. The present derivation is 
therefore equivalent to standard analyses. 

For our computation, we take the first cut imaginary part given by (||). 
The real part is found through an unsubstracted dispersion relation with 
77 = s/M 2 



sec 



0, 0, El) : 



Ti(M 2 ,x) = ReJi(M 2 ,x) 



M 3 x 



C<p 2 7 
M 2 



TT 



X 



ds 



4m2 (s - M 2 ) Z 



In 



1 + a/1 - Am 2 /s 
1 - a/1 - 4m 2 Js 



2 M 
— arctan — 
7r 27 



(10) 



where the first equality holds since M 2 < Am 2 . Remark that for x — ► 1, the 
second cuts vanish (see 6) and ([ji]) reproduce the result obtained in [11] for 
the parapositronium two-photon decay. 

Inserting the result for X\ into the expression of the differential rate with 

C 2 



M/m 2 and \<f) \ 

dT (p-Dm — » 7e + e 
dx 



a 3 m 3 1 



lir, we get: 



am 
6tt 



( Am 2 \ 


2 M 


\M 2 ) 


— arctan — 


7r 27 



p(x,a) 
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Note that, independently of the value for 7, the spectrum is always linear 
in x due to the 1/x in ©. This is an incorrect soft photon behaviour since 
Low's theorem requires a cubic spectrum instead [12] (when the amplitude 
behaves as x, the rate ~ x 2 and an additional x comes from phase-space) . 
By expanding the differential rate around 7 = 0, we get 



dT (p-Dm —* 7e + e ) 



am 



a 1 2 
1 - - + -a 2 
it 4 



dx 6ir 
As lowest order, we recover the standard result: 



dT (p-Dm — > 7e + e ) a 6 m M p (x, a) 



96vr ^ 64 ^ { ' x 2 



dx 



6?r 
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Note that the order a corrections disappear when using the form factor Fg 
(see [11]). For completeness, the total rate is simply : 



T (p-Dm — > 7e + e ) 



a 6 m 

67T 



F(a) 



( 4m 2 \ 


2 M 


\M 2 ) 


— arctan — 


7T 27 



with 



F(a) 



yi - a (a -4) + 2 log 



1 + yT 
1-vT 



3.3 The oblique cut contribution to the rate 

We have just seen that the vertical cuts suffice to reproduce the lowest order 
decay rate. It is interesting to investigate how the oblique cuts can restore 
the analytical behaviour of the spectrum without affecting this lowest order 
evaluation. 

The imaginary part corresponding to the oblique cuts is, from (6), 



ImX2 (s, x) 



s 16ttx 



dq 



87T7 



4m 2 



4m 



s 16ttx 

with g min , g max given in (7) and 



87T7 [H (g max ) - H (<?min)] s 



1 



H(q) 



2^ 



M(M 2 



In 



M + y/s + 2q 



+ 



M 2 



hi 



{^S + 2q) 2 -M 2 



We should now use an unsubstracted dispersion relation to get X2 (M 2 ,x). 
This integral is quite complicated and not very interesting for the present 
purpose. Instead, we revert to numerical evaluation of the dispersion integral 
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for T2 (M 2 ,x) = ReT2 (M 2 ,x) as a function of x, and compare it to (|T 
In the figure 3, we plot the vertical cuts contribution (dashed line) 

M 3 ^ /„,2 \ 2 M 
2i (M , x) = — arctan 



X ■ 



C4>o X ' vr 2 7 
and the complete result (solid line) 

M 3 _ , . x M 3 



X(M 2 ,x) = x- — - [Ti(M 2 ,x) + J 2 (M 2 ,x)] 



2 M M7 / ,+0 ° ds r „,, , „,, 
-arctan — + — / — ^ W g max (x ) - H (ft, 



for m = 1,M = 1.999 (i.e. 7 w 0.03) and for m = 1,M = 1.99994 (7 w 
0.0077). For x — > 0, one can verify that x -2" (M 2 , x) — > 0. The figures show 
that near zero the two cuts interfere destructively in order to maintain a 
correct analytical behaviour for the whole amplitude. Away from x = 0, the 
oblique cut contributions are strongly suppressed relatively to the vertical 
one, and this suppression increases as 7 decreases (I2 when 7 — ► 0). 
As can be seen on the graph, it is typically for x < 7 that the oblique cut 
contributes. Therefore, we can summarize by giving a simple representation 
of the different contributions to the amplitude. From the figure 3, one 
can see that the behaviours of the vertical cuts, the oblique cuts and their 
combination are quite precisely modelled as 

1 7/M f 1 

2i ~ — , Z2 ~ : 7— r =^ 2 = T\ + ±2 



x(x + 7/M) \x + j/M 

As a consequence, the spectrum behaves as 

dF(p-Dm^-fe + e-) 2 / % z ( 1 V f s 
- ■ |2] p (x, o) ~ x — — ~ x ' 



dx \x + j/M J \x + j/M 

(12) 

i.e. a linear spectrum when 7 — > (11), and a x 3 spectrum when x is small. 
The effect of 7 7^ is therefore to soften the photon spectrum and slightly 
reduce the total width. This behaviour is exactly what we postulated in a 
previous work [12|. 

In conclusion, the oblique cuts have a small contribution to the decay 
rate comparatively to the vertical cuts. However, their presence is essen- 
tial to guarantee the analytical properties of the amplitude expected from 
Low's theorem. Further, when tackling O (a 2 ) corrections, one must include 
contributions from the oblique cuts. 

4 Orthopositronium Decay to three Photons 

We now turn to the interesting decay o-Ps — ► 777. We start as usual with 
the loop model amplitude. Dispersion techniques express that amplitude 
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into separate contributions, arising from different cuts in figure 4. Since only 
off-shell intermediate leptons appear in the loop, the complete amplitude has 
a correct soft photon behaviour. This implies that one must consider all the 
cuts to preserve the analytical properties of the amplitude. This is one 
conclusion of the previous section. 

Let us give the loop model amplitude 

M^p (o-Ps -» 777) = e 3 e a (P) [ -^- A F B Tr J 7 a =i T^p r 

J (2vr) 4 B I 4-\P-m 4+k 



where the vertex is given by six amplitudes grouped as 

Y^p — y^ u p _|_ Y^p + T pl/p 

with 

Tt up = Y - >yP - ^ 

1 4-\P+V2-m 4+\P-Yi-m 

+ ryP _ ryP 1 ~» 

i-lP+Vx-m 1 4+\P-y 2 -m 1 

rf p = Y — 1 — t m — , — 7 P 

1 1 

+V V V 

i-lP+y^-m 1 4+\p-y 2 - m ] 

rf p = 7 M — i-^ i v — i-^ Y 

4- i p+ y 3 - m J 4+ \ p- y 1 - m 

and P = l\ + I2 + ^3- The amplitudes in T pup are shown on figure 4. Using 
charge-conjugation, we can show that the two drawn amplitudes are equal. 
The proof of this is in close analogy with the demonstration of Furry's 
theorem, and requires a change of integration variable q — > —q. This change 
is allowed for the form factor (see (||), (||)). Therefore, one can forget one 
term in each Y pvp and multiply the other by 2. In the following, to keep the 
combinational as clear as possible, we will continue to consider all the six 
diagrams. 

Let us now review the contributions aarising from the possible cuts in 
figure 4. 

4.1 The vertical cut reproduces standard approach results 

First we have the six vertical cuts shown on figure 5. As demonstrated in the 
paper pljl, those vertical cuts reproduce the known result: their combination 
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is what is calculated with (Q), i.e. at lowest order |3j 

, 2(tt 2 -9) fi 
T (o-Ps — > 777) = a m 

97T 

This result is obtained in the static limit, i.e. with Fg cx 5^ (q). 

If this limit is not taken, one has to go through the integration of the 
convolution- type amplitude (|3]). As explained in pi]] , the order of the cor- 
rections that will arise from that integration depends on the explicit form 
for Fb- Let us recall that for the parapositronium decay to two photons, we 
found that the form factors (||) led to O (a) corrections and higher, while 
(||) gave corrections starting at O (a 2 ) . 

Finally, it is important to note that the six vertical cuts are separately 
gauge invariant since the scattering amplitude for on-shell e + e~ — > 777 is. 

4.2 The oblique cuts and the structure dependent contribu- 
tions 

The oblique cuts are depicted in figure 6 for the photon 1 on the positron- 
ium side. Similar cuts for the photon 2 or 3 singled out are easily drawn. 
The oblique cuts are not gauge invariant, contrary to the vertical ones. 
To see this, it suffices to note that each of these oblique cuts is, from the 
optical theorem, a product of a scattering amplitude e + e~ — ► 77 times a 
bremsstrahlung amplitude o-Ps — > e + e~"f. Gauge invariance is broken by 
the bremsstrahlung amplitude, because the form factor is evaluated at dif- 
ferent momenta for different cuts. To visualize the situation, let us change 
the momentum parametrization and draw amplitudes for o-Ps — > e + e~7 (li) 
as shown in figure 7. Obviously, when contracted by the two amplitudes 
fail to cancel each other, due to the different momentum dependences of Fb 



M Breml = eF B (pi + h,p 2 ) \u(pi) 7 ^ i+ ^ _ m l a v fa) j e a (P) 

+eF B (pi,pa + h)\u(pi) 7" , \, 7 M « O2) 1 e Q (P) 

= ieF B ( Pl + h,p 2 ) U(pi) ^Ipji/ 1 ^ (Pa)} e « ( p ) 

+ie 2 F B (pi,P2 + h) («(pi) 7" ~y~f 17 ^ (P2)\ e a (P) 
{ 2p2 ■ h J 

and when contracted by : 

K„M Breml = ie [F B (pi + h,p 2 ) - F B ( P i,p 2 + h)\ {u (pi) ~/ a v (p 2 )} e a (P) 

Therefore, one should add to these diagrams the structure dependent am- 
plitude 7W^ D1 (figure 7). The new direct coupling between o-Ps,e + ,e~ and 
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7 must be such that the combination of the bremsstrahlung graphs and this 
structure graph is gauge invariant. For the paradimuonium, we have not 
considered such structure terms because the amplitude was gauge invari- 
ant throughout thanks to the factorized tensor structure (see (|j)). This is 
a peculiar feature of the two-photon decay of pseudoscalar positronium or 
dimuonium state. 

Let us now turn to the soft photon behaviour. The structure terms 
do not alter the cancellation among cuts in the soft photon limit (see the 
corresponding discussion for K$ decay in the appendix). Indeed, Low's 
theorem applied to the bremsstrahlung plus structure dependent amplitudes 
gives the following expansion around h = (see @, @): 



J ^ l Brem ^ Jvl Struct 



iF B (P1,P2) \ u(Pi 



v(p 2 )\e a (P)+0(h) 



2# + 7"fl Q r M+Vrf 

2p 1 ■ h 2p 2 ■ h 

i.e. terms of order cancel because -M^B rem + M Struct ls S au S e invariant. 
Importantly, note that the form factor is now evaluated at (p\ = m 2 ,p\ = to 2 ) , 
i.e. at the same point as for the vertical cuts. Since for a constant form fac- 
tor the combination of all the cuts necessarily behaves correctly in the soft 
photon limit, and since a momentum dependence in the form factor only 
modifies O terms, the conclusion follows. 

Remark that the appearance of these new structure contributions could 
have been guessed from the start, since the loop model amplitude (^) fails 
to be gauge invariant due to momentum dependences in the form factor. 
More precisely, the Ward identity l\^M^ vp (o-Ps — ► 777) = is verified 
only if the form factor allows for linear changes of the integration variable 
like q — * q + l\. A general form factor will not allow such shifts and one 
must supplement the loop model with new structure dependent amplitudes, 



as noted in [12]. The dispersion technique we followed here is interesting 
since we get a more precise information on their origin, and we are able to 
constraint them through Low's theorem. 

The same discussion can be made for the photons 2 and 3. One ends 
up with twelve oblique cuts, to which six structure dependent amplitudes 
must be added. All these new contributions should become important at 
the order 7 2 ~ a 2 /4, so that one can really question the completeness of the 
results given in the literature H. 



5 Conclusions 

We have shown that the standard approaches used to calculate positronium 
decay rates cannot be used to fully evaluate O (a 2 ) corrections. In the 
case of orthopositronium, many new contributions arise at that order from 
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processes where the electron emerging from the bound state decay is off- 
shell. Those contributions are essential to preserve the basic property of 
analyticity of the positronium decay amplitude, i.e. its vanishing in the soft 
photon limit. Further, gauge invariance requires some interplay between 
the annihilation process and the bound state dynamics through structure 
dependent amplitudes. 

All these new contributions vanish in the static limit, i.e. when the form 
factor is replaced by a delta function 6^ (q). Therefore, one can view these 
as non-perturbative effects arising from the binding of the constituents in the 
bound state. Indeed, the size of the corrections is fixed by the appearing 
in the bound state wavefunction. Concerning this wavefunction, forms where 
the dependence on the energy qo is replaced by a 5 (go) should not be used 
because it is specifically for non-zero value of the energy that the additional 
diagrams contribute. This raises the question of the form for Fb to be used. 



Let us recall that in [11], we advoquate the use of the improved form factor 
Fg instead of the usual Schrddinger momentum wavefunction. 

In conclusion, it should now be clear that O (a 2 ) corrections as given 
in the literature are rather incomplete. A lot of work is still needed before 
a definite theoretical prediction up to that order can be compared to ex- 
periment. Even if the theoretical basis seems settled, the orthopositronium 
life-time puzzle remains an open question. 
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6 Appendix : Kaon Radiative Decay Kg — > r ye + e 

The amplitude for Kg — > 7e + e _ is computed at lowest order in a pion loop 
model. The pion is treated as a point-like charged particle, which allows 
simple scalar QED treatment. Special attention is paid to the low energy 
behaviour of the amplitude. 



6.1 Decay Amplitude and Loop Integration 

The amplitude for kaon decay into / je + e~ is given by 

M(Kl^^e-) = -^M (K° - vr+vr-) £ * (k) {l1 {p) ]f {P ' )} 

x / -^ M>1V (* + *~ - rr) 

with t = p + p' = P — k. The 7r + 7r _ — ► 77 amplitude arises from the one- 
photon (figure 8a) and two-photon (the so-called seagull graph, figure 8b) 
coupling amplitudes as 

(2q - kf (2q + t) u - g» v (q 2 - m 2 ) 



(tt+tT -► 77) = 



({q + 1) 2 ~ m2 ^j ~ k) 2 ~~ m2 ) (i 2 ~ m2 ) 



with m the pion mass. The amplitude M (K° s — > 7T + 7r ), taken as a con- 
stant, has been factored out. The integration is done using dimensional 
regularization to preserve gauge invariance and get a finite result. We ob- 
tain 



I -^M» u (tt+tt- 77) 

y (27r) d 



" \*%L [g^ (k-t)- t»k v \ + idlzMg^f 
Jo Jo L 



(47r) 2 Jo Jo 

The denominator function is A = m 2 {1 — 4{a — b) xy + Aby (y — 1)) with 
the definitions a = M 2 /4m 2 ,b = t 2 /4m 2 and M the kaon mass. In the 
kaon rest-frame, b = M 2 (1 — oj) /4m 2 with oj the reduced photon energy 
2k°/M. When integrating over Feynman parameters, the last term vanishes. 
Therefore, the gauge invariance of the amplitude becomes manifest: 

— 2e 3 1 
M (ffg ^ 7e+e -) = (^tt+tt-) — 2 F (a, b) x 

^ (*) (P) 7„« (*/) } ^(k-^-^ (M) 
i 7 (a, 6) is the Feynman parameter integral 



/•l rl-y 

F(a,b) = dy dx, 
Jo Jo 



1 — 4 (a — b) xy + 4by (y — 1) + 
The prescription ie gives the sign of the imaginary part. 
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6.2 Feynman Parameter Integration via Dispersion Relations 



To calculate F (a, b), we shall use dispersion relations (see for example 14], 
[15|). A direct integration is possible, but we shall gain insight into the 
dynamics of the process by using dispersion techniques. 



Absorptive Part Extraction 

A first integration over the Feynman parameters gives 



F (a, b) 



dy 



y-1 In K(b)- InK (a) b (y - 1) (InK (b) - InK (a)) 



Ay {a -by 



{a -by 



, (15) 

where In K (x) = In (1 + 4x (y — 1) y — ie). The imaginary part of F (a, b) 
comes from the values of y for which the argument of a logarithm is negative. 
There the logarithm imaginary part is —in. Integrating on the relevant y 
values, we find 



ImF(a,b) 



7r (in [-y/a + \/a — 1 



bJ^ 



-/r ( In 



2 (a - b) z 
Vb + Vb^T 



2 {a -b) z 

or, recalling the definition of a and b (see 



(a-1) 
■9(b-l) 



(16) 



ImF(a,6) 



M 4 



hi 



I + a/1 



4m 2 
A? 2 " 



4m 2 
A? 2 " 



27rm 2 1 — uj 
~M~ 2 



4m z 
TT 1 



9 (M 2 - Am 2 ) 



47rm 
M 4 w 2 



In 



I + a/1 



4m. 2 
A/ 2 (l-u;) 



l 



4m 2 
A/ 2 (1-uj) 



1 



2-irm 2 I — uj 
~M~ 2 ^2" 



4m 



(17) 

Let us consider the figure 9. Obviously, the first line in (17) corresponds 
to the vertical cuts (which contribute only if the decaying particle's mass 
M 2 is greater than (2m) 2 ) while the second line corresponds to the oblique 
cuts (which contribute only if the virtual photon energy t 2 is greater than 
(2m) 2 ). 

Further, one can see that in the soft photon limit uj — > (or b — > a), 
ImF (a,b) behaves as a constant. Indeed, applying L'Hospital's rule twice, 
we get 



ImF ( M 2 M 2 (1-uj) \ ^ 
yAm 2 ' 4m 2 J 



2m q 



M A y/l-4m 2 /M 2 



(18) 
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while individually the contribution of each cut is divergent as l/(a — b) 2 ~ 
1/u 2 . We will come back to this point later. 



Dispersive Part Integral 

In order to write down the unsubstracted dispersion integral, let us ex- 
press ImF (a, b) as a function of the available energy through a (s) = s/Am 2 
and 6(s) = s (1 — u>) /Am 2 . Then 

Re F (a (s Q ) , b (s )) = - / — ^- Im F (a (s) ,b(s)) (19) 

IT J S- S 

with s D < 4m 2 , such that we can omit the principal part. For such a 
kinematics F (a (s a ) , b (s G )) = ReF (a (s Q ) , 6 (s G )). In the next section we 
shall analytically continue F to the physical value s a = M 2 . The result of 
(19) is easily obtained in terms of the integrals 

1 dy n ri+v^l o ■ 2 P~ f 1 d v 



In 



o z/o - y 



2arcsin 2 W — , / -i/l — y = 2—2\/ y a — 1 arcsin 

w 2/o Jo Vo-V 



valid for y > 1 (see JlJ]). We can write F (a,b) with < a < 1 and 
< b < 1 as 



^ fa & ) = "TTT^T + 7 772 ( o (/(«)" / (*0) + % W - 5 (6)) ) (20) 



1 1 fl 

2(a"^6) + ( a -6) 2 \2 

in terms of / (x) = arcsin 2 y^x and g (x) = yj arcsin y/x. 
Analytic Continuation 

The results for a > 1 and b > 1, are obtained by analytic continuation of 
(20). The analytic continuation of / (a) and g (a) are 

/ (x) = arcsin 2 {\fx) = — ^ln (yfx + \J x — l) — 2* 7r J 

5 (x) = \J X arcsin \/x = \J X ^ln (y/x + \fx — l) — T^ 71 "^ 

The equalities hold for x G C, Imx > 0. The right hand sides define the 
analytical continuation for x > 1. Hence the result is conveniently expressed 
as © with 

/ arcsin 2 (^x) < x < 1 

/W 1 -(ln^+v^l)-!^) 2 x>l 1 j 



0*0 



arcsin (-^/i) < x < 1 

p (In (-v/x + y 7 ^ - l) - |wr) x > 1 



x-1 
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Had we chosen to analytically continue / and g in the lower half complex 
plane, their imaginary parts would have had the opposite sign. One can 
verify from (pll) that the upper half complex plane analytic continuation 



reproduces ImF (a, b) as given in (16). This result corresponds to [ 20 1 where 
a sign mistake has to be corrected, and to pl|, p2]. 



6.3 Decay Width, Differential Width and Low's Theorem 

The differential decay width is given by dT (K s — > e + e~j) = Ylspin \M 1 2 
After a straightforward integration over the electron-positron phase-space, 
we get the differential rate in terms of the reduced photon energy u: 



dF (K s -> e+e- 7 ) /duj _ a 3 



F 



1 l-u 



"" [a e + 2(l-w)] 



T (K° s -» vr+vr-) 3^ 3 o 2 ^!^^ \ 1 to v n (1 

with a e = Am 2 /M 2 and = 4m 2 /M 2 = 1/a. The bounds on the u 
integration are w m i n = and u; max = 1 — a e . A numerical integration of the 
differential rate gives the prediction for the rate Kg — > e + e~7 relatively to 
K s -> tt+tt- 



T (K s -» e+e-7) 



~ 4.70 x 10 -8 (22) 



77 w ~ " r [k° s 

where the real and imaginary part of F ( — , — — ) contribute for 1.26 x 10 -8 

and 3.43 x 10 -8 , respectively. 

Let us turn to the soft photon behaviour of the differential width. We 

have to analyze the loop integral function F ■ This function tends 

to a constant for very low u: 

1 1 — aA a w a 2 



(note that for a w < 1, the imaginary part of (23) is the same as in (|18|)). 
Let us emphasize the strong cancellation between the two cuts in the soft 
photon limit. Individually, each cut diverges as to 2 , but their combination is 
convergent. Since the rest of the amplitude tends to zero like u> (due to the 
[g^ u (k ■ t) — t^k v ] factor), Low's theorem is verified: the amplitude behaves 
as to near u = 0. The resulting spectrum is thus in lu 5 (lo 2 from the squared 
amplitude and a uj from phase space). In other words, if we had forgotten 
one cut, the resulting spectrum behaviour would have been divergent as 1/uj 
near zero instead of vanishing like to 3 . In fact, this a; 3 spectrum is exactly 
the same as in the ir° — ► e + e _ 7 differential decay rate. 
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Comparison with the Two-Photon Decay Mode 



From the amplitude (|14"|), we readily obtain the amplitude for — > 77 
by removing the electron current, the photon propagator and by taking the 
limit u — ► 1: 



M(K° S 



77) 



-2aM(K° s ^tt+tt~) # 



4-7T 



{k)e* u (k') [g^(k-k') -k'^k v ]F 



withF(i )=K{/(i)-i}i-e. 



?7T 



Am 



M 2 \ M 2 



hi 



1+A/1- 



i-a/i- 



4 m 2 




7T 1 , 2 

1 In 

2 M 2 M 2 



which corresponds to the result given in 
rate is : 



1 21], [23]. The relative decay 



77 



a 



1 



vr 2 a 2 VT 



F I — ,0 



2.94 x 10 -6 



while the experimental value for this ratio is (3.5 ± 1.3) x 10 6 [24] 
We can now compare the two electromagnetic modes : 



R 



e T e 7) 



77 



0.016 



r [k° s -> 77) 

in accordance with [21] and Jl9| ]. This ratio i? 77 is similar to [24] 

T (vr — ► e + e~7 N 



r (-7T — ► 77) 



0.012 



The small difference is due to the phase-space factor. The similarity between 
a constant coupling model like for 7r° — ► e + e _ 7 and the present loop model 
can be understood from the behaviour of the photon energy spectrum. In- 
deed, the pure phase-space spectrum is very strongly peaked at high to, and 
therefore the decay rate is quite insensible to the detail of the F Q^, 

function. Specifically, if we replace F 



1 l-w 



by its value for u> = 1 (|24|) 



we find the ratio iL 



~ 4.671 x 10 , very close to 



6.4 Soft Photon Behaviour with a non-constant Form Factor 

We now turn to the case of a non-constant amplitude A4 (Kg — > 7r + 7r~) = 
T ((q + k) 2 , (q — t) 2 ^j . We analyze the process from the point of view of 
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Low's theorem. The constraint of gauge invariance will be seen to require 
some new structure dependent contributions. We will prove that those new 
contributions do not affect the interference between absorptive parts ob- 
served in previous sections, responsible for the vanishing of the amplitude 
in the soft photon limit. 

Let us concentrate on the absorptive part of the amplitude for K$ — ► 
7e + e~. Using the optical theorem (see figure 2), we write 

21m M^(K% -^je + e~) = V 1 + V» (25) 

with V M the vertical cut contribution and V* 1 the diagonal one : 

V = Jd*„ n (27r) 4 5^(P- Pl - P2 )M(K s ^7r + (p 1 )7T-(p 2 )) 

xM^ 1 (tt + (jpi) 7r~ (p 2 ) — > 77* — ► 7 e + e~) 
V» = J d*^ (2vr) 4 5^ (P-P1-P2- k) {K° s - vr+ ( Pl ) ir~ (p 2 ) 7 (fc)) 

X.M (7r + (pi) 7r~ (p 2 ) -> 7* -»■ e + e~) 

The 7T7T — > 7e + e _ amplitude is gauge invariant, as can be easily verified, 
and so is V M , independently of the form of T (pi,P 2 )- Therefore V^ 1 has to 
be gauge invariant too. The two bremsstrahlung contributions to T>^ are 

Ml B {K° s ^ir+( Pl )iT- (p 2 )j(k)) = 

- P» 2 



T ((pi + kf ,p 2 ) - ^ (p?, (P2 + fe) 2 ) 



pi • k V / p 2 ■ k 

Trivially, gauge invariance is verified if the form factor is constant. Problems 
arise if F (pi + k,p 2 ) ^= T (p±,p 2 + k), since when contracted by fc M the two 
terms fail to cancel. To maintain gauge invariance, one must supplement 
the amplitude Kg — > 7r + 7r~7 with a structure dependent term, for example 
M» D (K s ^ir+(p l )Tr~(p 2 ) 1 (k)) = 

> (K + *) 2 ) " ^ ((Pi + kf ,pl)} {P \~ P2) \ (26) 
LV / V / J pi ■ k — P2 ■ k 

such that fc".M£ B + k»M^ D = 0. 

Let us now analyze the soft photon behaviour of this process. The form 
factor is expanded around k = as 

^((Pi + fe) 2 ,p|) = ^(p 2 ,p2)+2p 1 -A;x^ / + 0(A; 2 ) 
^ (pi, (P2 + fc) 2 ) = ^(p 2 ,p^)+2p 2 -A;x^ , + 0(fc 2 ) 



since we expect T 1 = (Pi^P 2 ) = (Pi'P 2 ) when evaluated on-shell, i.e. 



at p 2 = p 2 . = "i 2 - The amplitude is now 



^ ib+sd = f{pIpI) 



Pl P 2 



Pi - k p 2 - k 



0{k) (27) 
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i.e. the terms of order (k)° disappeared. This result is valid in full generality, 
it does not depend on the present structure dependent term (26) chosen for 
illustration (see O, pi). 



Now let us analyze the soft photon limit for the imaginary part (25). We 
have obtained for the two cuts the structures 



V = T{pl P l) 
2?" = F{pl,p 2 2 ) 



T 



k 



+ O (k) 



Except for the O (k) terms in T>^, this is exactly the soft photon expansion 
one would get for a constant form factor. Hence we expect from the analyzes 
of previous sections that the combination will vanish when k — * 0. Note that 
it is essential that no constant term (due to the momentum dependence 
of T) appears in ■Mj B+SD . This was guaranteed by the gauge invariance 
of ■M'ib+sd^ w hich in turn necessitates the contribution of some structure 
dependent amplitude. 
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Paradimuonium Decay to ye e" : Figures 




Figure 2 : The vertical and oblique cuts contributing to the amplitude imaginary part 




Figure 3 : The vertical cuts (dashed line) and vertical plus oblique cuts (continuous line) 
contributions, for y=0.03 and y=0.0077 as a function of the photon energy x 



Orthopositronium Decay to three Photons : Figures 




Figure 4 : The loop model diagrams r; 




+ 5 photon permutations 



Figure 5 : The vertical cuts, known to reproduce standard results for the o-Ps decay rate 




Figure 6 : The oblique cuts. Some new 0[a) contributions to the o-Ps decay rate 



Orthopositronium Decay (contd.) : Figures 




Figure 7 : The Bremsstrahlung and strupture amplitudes, contributing to the oblique cuts 



Kaon Decay to ye e" : Figures 




Figure 9 



The vertical and oblique cuts 



contributing to lmF(a,b) 



